No Title  |POLVER05_3
d(Ca)/d(t) = -r1 #tcpo
d(Cb)/d(t) = -r1 #h2o2
d(Cc)/d(t) = (2*r2) #PhOH
d(Cd)/d(t) = (r1 - r2) #peroxyacid
d(Ce)/d(t) = (r2 + r6 - r3) ## 1,2-dioxetanedione
d(Ch)/d(t) = -2 * r6 + 2 * r3 #CO2
d(Cf)/d(t) = r5 + r4 + r6 - r3 #dye
d(Ck)/d(t) = r3 - r4 - r5 - r6 #complex
r1 = k1 * Ca * Cb #
r2 = k2 * Cd #
r3 = k3 * Cf * Ce #
r4 = k4 * Ck #
r5 = k5 * Ck #
r6 = k6 * Ck * (Ch) ^ 2 #
rEMIT = r4 * QY #light as a function of time
QY = 200000 #quantum yield (overall height)
k1 = 1.485 * 10 ^ 7 #
k2 = 1.485*10^7 #
k3 = 891 #
k4 = 0.05 #
k5 = 0.111 #
k6 = 1.782*10^12 #
t(0)=0
Ca(0)=0.7
Cb(0)=0.8
Cc(0)=0
Cd(0)=0
Ce(0)=0
Ch(0)=0
Cf(0)=0.0005
Ck(0)=0
t(f)=3600